Can one count the shape of a drum? 
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Sequences of nodal counts store information on the geometry (metric) of the domain where the 
wave equation is considered. To demonstrate this statement, we consider the eigenf unctions of the 
Laplace-Beltrami operator on surfaces of revolution. Arranging the wave functions by increasing 
values of the eigenvalues, and counting the number of their nodal domains, we obtain the nodal 
sequence whose properties we study. This sequence is expressed as a trace formula, which consists 
of a smooth (Weyl-like) part which depends on global geometrical parameters, and a fluctuating part 
which involves the classical periodic orbits on the torus and their actions (lengths). The geometrical 
content of the nodal sequence is thus explicitly revealed. 
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Eigenfunctions of the Schrodinger and other wave 
equations can be characterized by the number of their 
nodal domains - a nodal domain being a maximally con- 
nected region where the eigenfunction has a constant 
sign. The intimate connection between the spectra of 
wave equations and the corresponding sequences of nodal 
counts is well known and frequently used in various 
branches of physics and mathematics. Sturm's oscilla- 
tion theorem states that in one dimension the n-th eigen- 
function has exactly n nodal domains. In higher dimen- 
sions Courant proved that the number of nodal domains 
v n of the n-th eigenfunction cannot exceed n |]|. Re- 
cently, it was shown that the fluctuations in the nodal 
sequence {v n }^\ display universal features which distin- 
guish clearly between integrable (separable) and chaotic 
systems Their study also leads to surprising connec- 
tions with percolation theory Q. Moreover, the nodal 
sequences of several isospectral (yet not isometric) do- 
mains were recently shown to differ in a substantial way 
0, . The later observations suggest that the nodal se- 
quence stores information about the domain geometry, 
and this information is not equivalent to the one stored 
in the spectrum. Here, we provide further evidence by 
deriving an asymptotic trace formula for the nodal count- 
ing function 
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K > ; [•] = the integer part. (1) 



The trace formula (see (J7J) and (|19|l ) shows the explicit 
dependence of the nodal sequence on the geometry of 
the surface in both the smooth (Weyl-like) and the fluc- 
tuating parts. Thus, the nodal trace formula is similar 
in structure to the corresponding spectral trace formula 
0,0,0. Kac's famous question "Can one hear the shape 
of a drum?" was triggered by the study of the progeni- 
tors of the spectral trace formulae Q . The trace formula 
for the nodal counts leads us to the title of this letter in 



which "count" replaces "hear" . We will consider here two 
particular classes of systems, namely, the wave equation 
on convex smooth surfaces of revolution and on simple 
two-dimensional tori. Generalization to other Rieman- 
nian manifolds in two or more dimensions are possible, 
provided the wave equation is separable. 

The nodal counting function Q is well defined if 
the spectrum is free of degeneracies, E n > E n —i- In 
case of degeneracies we represent the wave functions 
in the unique (real) basis in which the wave functions 
appear in product form. This, however, does not suf- 
fice to set a unique order within the degenerate states, 
which consequently introduces ambiguities in the nodal 
sequence. To circumvent this problem, we modify the 
definition of the nodal counting function: First define 
c(-E-) = J2^=i v n < d{E — E n ). This function is based 
on information obtained from the nodal sequence and 
the spectrum. To eliminate the dependence on the lat- 
ter, we use the (e-smoothed) spectral counting function 
Af e (E) = E£Li(i + h arctan E =^), which for finite e is 
monotonic and can be inverted. We define E e (K) as the 
solution of Af e (E) — K and the modified nodal counting 
function is 
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If there are no degeneracies, c(K) is equivalent to 
up to a shift K — > K — ^. A g-times degenerate eigen- 
value E n = E n+ i = ■ ■ ■ = E n+g ^i contributes a single 
step function 6 (if — (n — 1 + §)) J2s=i u n+s-i where 
the nodal counting function increases by the sum of the 
nodal counts within the degeneracy class. We will derive 
a trace formula for this modified nodal counting function 
below (and omit the 'modified' in the sequel). 

We start with the simpler case of a 2-dim torus repre- 
sented as a rectangle with side lengths a and b and peri- 
odic boundary conditions. The eigenvalues take the val- 



ues E„ 



(27T) 



where m,n S Z. The cor- 



responding wavefunctions for m,n > are ipn,m(x,y) 



2 



cos(27rna;/a) cos(2irmy / b) and the cosine is replaced by a 
sine for negative m or n. The number of nodal domains 
in the wavefunction ipn.m is 



(2|n| + <5 n , )(2|m| + <W)) 



(3) 



The aspect ratio r — a/b is the only free parameter in this 
context because the number of nodal domains is invariant 
to re-scaling of the lengths. 

Using Poisson's summation formula J^"L n f(n) = 

Ew=-oo Ini-i/2 dn f(n)e 2 ™ Nn the leading asymptotic 
trace formula for the spectral counting function Af(E) = 
E m ,„ ®( E - En,m) for a torus, is, 

N{E) = AE (4) 



+ 



Here, „4 = ab/(Aw), and the sum is over the winding 
numbers r = (N,M) £ Z 2 \(0,0) (in the sequel every 
sum over r will not include (0, 0) unless stated other- 
wise) . L r = J (Na) 2 + (Mb) 2 is the length of a periodic 
geodesic (periodic orbit) with r = (N,M). 

Our goal is to derive a similar trace formula for the 
leading asymptotic behavior of the nodal counting func- 
tion. Again, using Poisson re-summation and the saddle 
point approximation we get for c(E) 

2A 2 
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To express the counting function as a function of the in- 
dex K, we formally invert the spectral counting function 
to order 0(K°) 
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where l T = L r /VA is the re-scaled (dimensionless) length 
of a periodic orbit. This formal inversion needs a proper 
justification which makes use of the fact that we actually 
invert the smooth and monotonic Af e (E). However, a 
detailed discussion of this point goes beyond the scope of 
the present note. The numerical tests which we provide 
here, supports the validity of this formal manipulation. 
Replacing E by E(K) in © and keeping only the leading 
order terms, we get the nodal counting function, which 
we write as a sum c(K) = c(K) + c osc (K) of a smooth 
part c and an oscillatory part c osc : 

c(K) = ^-K 2 + 0{K) , 

c osc (K) = ifi^a r sin(/ r V^- j) + 0(K) (7) 
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While the smooth part is independent of the geometry of 
the torus, the oscillating part depends explicitly on the 
aspect ratio r and can distinguish between different ge- 
ometries. The main difficulty in computing higher order 
corrections to the leading behavior of the nodal count- 
ing function is that products of sums over periodic orbits 
appear already in the terms of order K. 

Turning now to more general surfaces we consider ana- 
lytic, convex surfaces of revolution M. which are created 
by the rotation of the line y = f(x),x £ I = [— 1, 1] about 
the x axis. To get a smooth surface, f(x) in the vicinity 
of x = ±1, should behave as f 2 (x) w a±(l =F x), with a± 
positive constants. Convexity is achieved by requiring 
the second derivative of f(x) to be strictly negative, and 
f'(x m ax) = 0, x max € I, where / reaches the value / max . 
We consider the wave equation — Aip(x,9) = Eip(x,9) 
and the Laplace-Bcltrami operator for a surface of revo- 
lution is 
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f{x)o(x) dx a(x) dx f(x) 2 dd 2 



(8) 



Here, a(x) = + f'(x) 2 and 8 is the azimuthal an- 
gle. The domain of A are the doubly differentiable, 2tt 
periodic in 9 and non singular functions on [/ x S ]. 
Under these conditions, A is self adjoint and its spec- 
trum is discrete. A is separable and the general solution 
can be written as a product ^(cc, 0) = exp(im#) (j)m(x) 
where rn e Z. For any m, (jHJ) reduces to an ODE of the 
Sturm-Liouville type, with eigenvalues E m , n (doubly de- 
generate when m 7^ 0) and eigenfunctions (j) n ,m(x) with 
n = 0, 1, 2, . . . nodes. The eigenfunctions corresponding 
to the eigenvalue E n ^ m can be written as linear combi- 
nations of cos(m9)(j) nim (x) and sin(m9)(j) nim (x) . To be 
definite, we chose these two functions as the basis for the 
discussion, and associate the former with positive values 
of m and the later with the negative values of m. The 
nodal pattern is that of a checkerboard typical to sepa- 
rable systems and contains 



v n ,m — (n + l)(2|m| + <y m ,o) 



(9) 



nodal domains. The semi-classical spectrum is con- 
structed by using the Bohr-Sommerfeld approximation 

0, 

E s £ m = H{n+-,m) + G{l) ,n e N, m e Z . (10) 

where H (n, m) is the classical Hamiltonian defined in 
terms of the action variables m and n, where m is the 
momentum conjugate to the angle 9 and n is 

n(E; m) = J- (p p x (E 7 x) dx = - [ p x {E, x)dx 

27r J n j x _ ( 1:L ) 

p x {E, x) =^/(Ef 2 (x) - m 2 )(l + f'(x) 2 )/f(x) . 

x± are the classical turning points Ef 2 (x±) — m 2 = 0, 
with x_ < 

•Emax — •£+ • The h&miltonictii is obtained 
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by expressing E as a function of n, to using the implicit 
expression . H(n,m) is a homogenous function of 
order 2: H(Xn, Am) = X 2 H(n,m). It suffices therefore 
to study the function n(m) = n(E = 1, to) which defines 
a line T in the (n, m) plain. n(m) is defined for \m\ < 
/max and is even in to, n(m) — n(—m). The function 
n(m) is monotonically decreasing from its maximal value 
n(0) to n(m — f max ) = 0. All relevant information on 
the geodesies on the surface can be derived from n(m). 
Periodic geodesies appear if the angular velocities u> n — 

are rationally related. Since 



dH(n.m) 

dn 
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dH(n.m) 



fois is equivalent to the condition 



M + N 



dn(m) 
dm 
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for M,N ^ 0. The integers r = (M, N) G Z 2 \(0,0) are 
the winding numbers in the 9 and x directions. The clas- 
sical motion is considerably simplified if the twist con- 
dition 8] n"(m) = d for < to < / max is 
fulfilled. This excludes, for example, the sphere but in- 
cludes all mild deformations of an ellipsoid of revolution. 
We will assume the twist condition for the rest of this 
letter. It guarantees that there is a unique solution to 
(|12fl which we will call m r . Note, that n'(m) has a finite 
range CI and a solution only exists if —M/N £ CI. The 
cases N = 0, M ^ or M ^ 0, TV = are not described 
by solutions of l)12|l . They describe a pure rotation in the 
9 direction at constant x — x max where Tio,±|M — ±/max 
(N = 0) or a a periodic motion through the two poles at 
fixed angle 9 mod it (M = 0) such that mijv|,o = 0- The 
length of a periodic geodesic is given by 



2-7T \Nn(m r ) + Mm T 



(13) 



Returning to the spectrum, we note that the leading 
terms in the trace formula for the spectral counting func- 
tion N(E) = J2m n ©(-S - En,m) can be obtained by using 
(110(1 and Poisson's summation formula || . 



M(E) = AE + EiJ2^( E ) 



where 
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n(m)dm = ||.M||/47r 



(14) 



(15) 



and 1 1 M\ | is the area of the surface. The oscillating parts 
contain integrals cx /f" ax dm e 2mE ^ (Nn(m)+Mm) which 

can be calculated to leading order in E% using the sta- 
tionary phase approximation. The points of stationary 
phase are identified as the classically periodic tori l|12|) 
with to = TO r . This restricts the range of contributing r 
values to the classically accessible domain —M/N E O. 
Thus, 



2TT\N 3 n'^\^ 



0{E- 



(16) 



where n" = n"(m = m r ) and a = sign(n^) which is the 
same for all values of r. The contributions of the terms 
with cither N = or M = or with —M /N £ CI are of 
higher order in 1/E and will not be considered here. 

We are now ready to derive the asymptotic trace for- 
mula for 



c{K) = c(E(K)) = 
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The second equation was obtained by inverting K — 
Af(E) to the desired order using the trace formula (THf . 
We follow the same approach as for AT and expand 
the result in SE = E(K) - K/A such that c(K) = 
c(K/A) + c'{K/A)5E + 0{c"5E 2 ) which is consistent if 
we neglect all orders smaller than 0{K). The result can 
be expressed as a sum c(K) = c(K)+c osc (K) of a smooth 
part c and an oscillatory part, c osc , in complete analogy 
to 0. 
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Z (K) = K z y a r sin(l r K 
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(18) 



m p n q = — I dTOdn |TO| p n 9 

' E(rn,n)<l 

where l r is the re-scaled length of a periodic geodesic 
and a r is the amplitude contributed by the (classically 
allowed) r torus. For m r = or m r = ±/ max only one 
half of the stationary phase integral contributes and the 
amplitude a r has to be multiplied by 1/2. 

The approximations involved in the above calculation 
have been tested on an extensive numerical data base 
for two ellipsoids of revolution defined by the equation 
f(x) = — x 2 (R = 2 in data set a), and R = 1/2 in 
data set b)) and for two different tori (r 2 = 2 in data set 
c), and t 2 = v2 in data set d)). The spectral interval 
used for the numerical tests included the first 10 5 eigen- 
values for the ellipsoids, and the first 4 x 10 6 eigenvalues 
for the tori. The numerically computed c(K) were fitted 
to a fourth order polynomial in n = \[~K and in all cases, 
the agreement of the two leading coefficients with the 
asymptotic theory was better than a percent. The oscil- 
lating part has been obtained numerically by subtracting 
the best polynomial fit from the exact c(K). The fluctu- 
ating parts of the trace formulae were tested in two ways. 
The integral of the squared oscillatory part 



R(K) 



K 



dK' c osc (K'f 



(19) 
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FIG. 1: Numerical checks of the fluctuating parts of trace formulae for the two ellipsoids ( a) and b), see text) and the two tori 
( c) and d), see text), left: Double logarithmic plot of the integrated squared fluctuations 121 H (arbitrary scale), the full line 
has slope 7/2. Right: Length spectra of the nodal counting functions 12111 . The full line is obtained from the trace formulae 
1191 and Q and the points represent the numerical data. 



was computed as a function of K, and compared with 
the theoretical expression which consists of a double sum 
over periodic geodesies. Simplifying this expression by 
considering only its diagonal part, one gets the estimate 



(20) 



used to resolve isospectrality. Thus, the geometrical in- 
formation stored in the nodal sequence is not equivalent 
to the one stored in the spectral sequence. This result 
together with the trace formula obtained here motivates 
further research of the nodal sequence as a tool in spec- 
tral analysis. 



which scales like K' 7 1 2 . This scaling has been tested and 
the results are shown in the left part of figure B Clearly, 
the expected power law is reached for sufficiently large 
values of the counting index K. A more stringent test 
of the trace formula is provided by computing the length 
spectrum, defined by the properly scaled Fourier trans- 
form of c osc (K) with respect to k = \f~K. 



S(l) = I 3 / 2 



d/t K 5/2 c n 



Z (K = n 2 )e 



-\-ikI 



(21) 

Gaussian windows centered at k = kq and with a width 
oc y/ui restricted the data used to be well within the semi- 
classical domain. The trace formula for the nodal counts 
predicts pronounced peaks at the scaled lengths I = l r of 
the periodic geodesies. The right frame in Figure B shows a 
remarkable agreement of the numerical data with the the- 
oretical predictions. This excellent agreement provides 
further support for the validity of the approximations 
which were used in the derivation of the two versions of 
the nodal counts trace formula. 

Recent studies 0, have shown that the nodal se- 
quences of isospectral domains are distinct, and can be 
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